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We report numerical estimates of correlation lengths in 2D Potts models from the asymptotic decay of the 
cluster-diameter distribution. Using this observable we are able to verify theoretical predictions for the correlation 
length in the disordered phase at the transition point for q = 10, 15, and 20 with an accuracy of about 1% — 2%. 
This is a considerable improvement over previous measurements using the standard (projected) two-point function. 



1. INTRODUCTION 

The two-dimensional g-state Potts model [1] is 
one of the most thoroughly studied models of sta- 
tistical physics. Since many quantities of interest 
are known exactly [2] it has repeatedly been used 
as a testing ground for more general ideas on, e.g., 
finite-size scaling (FSS) [3], the shape of energy or 
magnetization distributions [4] or partition func- 
tion zeros [5]. A few years ago also an analyt- 
ical formula for the correlation length £d(Pt) in 
the disordered phase at the first-order transition 
point p t for q > 5 has been added to the list of 
exactly known quantities [6,7]. For the correla- 
tion length £o{fit) in the ordered phase, however, 
only heuristic arguments were available until very 
recently [8]. 

This motivated us to perform a numerical study 
with the goal to clarify conflicting predictions for 
the ratio £ /£d at (3 t [9]. To test the employed 
numerical techniques, we started this project by 
first studying the disordered phase [10]. Using 
cluster algorithms for the update and standard 
(projected) two-point correlation functions for ex- 
tracting the correlation length we could confirm 
the analytical formula only with an accuracy of 
about 10% — 20%, depending on the value of q. 
The deviation is clearly of systematic nature and 
could be traced back to the unexpected impor- 
tance of higher order excitations. Even though 
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we worked on large lattices of linear size up to 
« 28£d and performed quite a large number of 
about 50 000 - 100 000 statistically independent 
measurements, fits with sufficiently many correc- 
tion terms turned out to be too unstable. 

Here we report a reinvestigation of this prob- 
lem by using a completely different estimator for 
the correlation length, the asymptotic decay of 
the cluster-diameter distribution as discussed in 
recent analytical work [8]. Numerically this ob- 
servable turns out to be very well suited in the 
disordered phase, allowing for the first time a con- 
firmation of the analytical formula for £d(Pt) with 
an accuracy of about 1% — 2%. 

2. SIMULATION 

We simulated the standard Potts model parti- 
tion function 

z = E e ~ pE ~> E = ~ E ; * = i, ■ ■ • . <*. (i) 

where i denote the lattice sites, (ij) are nearest- 
neighbor pairs, and 5 SiS . is the Kronecker delta 
symbol. All simulations were performed in the 
disordered phase at the infinite volume transition 
point pi = \n(l + y/q) for the cases q = 10, 15, and 
20. As in our previous studies [10] we minimized 
autocorrelation times by updating the spins with 
the single-cluster algorithm. To take advantage of 
translational invariance we always employed pe- 
riodic boundary conditions. By using very large 
lattices of size V = L x L with L = 300, 120, 
and 80 for q = 10, 15, and 20, respectively, we 
achieved that, starting from a disordered config- 
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uration, no tunnelings into the ordered phase oc- 
curred during the simulation. The code was im- 
plemented on a T3D parallel computer in a triv- 
ial way by running 64 time histories in parallel. 
In this way we generated a very high statistics 
of 600 000 T int , e sweeps for q = 10, 9 000 000 T int , e 
for q = 15, and 4 200 000 T int , e for q = 20, where 
Ti„t, e (~ 59, 18, and 25) is the integrated autocor- 
relation time of the energy. The error bars are 
estimated by jack-kniving over this statistically 
independent sample. 

In this work we mainly focussed on the cluster- 
diameter distribution G (x) which is defined 
[8] as the probability 



G (x) = fi( diam C tQ = x) 



(2) 



that the cluster C tQ connected to a lattice site 
io has a given diameter x, where the cluster di- 
ameter diam C; is the maximal extension in any 
of the D coordinate directions; see Fig. 1. Since 
we employ periodic boundary conditions, we have 
averag ed G (x) over all lattice sites io- In prac- 
tice this amounts to histograming the cluster di- 
ameter weighted by the cluster size. The theoret- 
ically expected asymptotic behaviour of G (x) 
is an exponential decay governed by the correla- 
tion length £d [8], 



G (x) = aexp(-a;/£d). 



(3) 



Taking the logarithm of G and performing lin- 
ear two-parameter fits it is then straightforward 
to extract 

For comparison we considered also in the new 
simulations the ky"^ = 2irn/L momentum projec- 
tions (i = (i x ,i y )), 



iky ^ (*y Jy ) 



of the spin-spin correlation function 

G(i,j) = (s titj -^) = ^zl{e(i,j)), 



(4) 



(5) 



where the last equality gives the improved clus- 
ter estimator with Q(i,j) = 1, if i and j belong 
to the same cluster, and 6 = otherwise. As 
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Figure 1 . Illustration of the definition of the clus- 
ter diameter diamCi = max{£ x ,£ y }. 



discussed previously [10], here non-linear four- 
parameter fits of the form 

gW (x) = a ch( L/2 ~ X ) + b ch(c L/2 ~ X ) (6) 
are necessary to extract 

3. RESULTS 

In all simulations we monitored the energy and 
magnetization to convince ourselves that the sys- 
tem never tunneled into the ordered phase. As a 
more quantitative measure we also recorded the 
maximal cluster diameter and computed energy 
and magnetization moments which can be com- 
pared with exact or series expansion results. As 
a result we are convinced that, despite the very 
long run times, our results for ^ can be identified 
with the pure phase correlation length. 

The data for G """(x) and g^°\x) are shown 
for q = 10 in the semi-log plot of Fig. 2. The con- 
tinuous lines are one- and three-parameter fits to 
the Ansatz (3) and (6), respectively, with ^ held 

fixed at its theoretical value 10.5595 We see 

that this yields an excellent description of the fall- 
off of g(°^ over more than four decades. Still, from 
a four-parameter fit over the same x range with 
£d as a free parameter we estimate an about 10% 
smaller value of ^ = 9.5(4). This confirms our 
earlier observation that four-parameter fits to 
systematically underestimate Compared with 
Ref. [10] the statistics of the present simulations 
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Figure 2. Semi-log plot of the cluster-diameter 
distribution G (x) and the projected correla- 
tion function g^ for q = 10 at [i t in the disor- 
dered phase. 



is higher by more than one order of magnitude. 
This allowed us to include larger x values in the 
fits and, as expected, improved the estimates of 
£d, in particular for q = 15 and 20. An attempt to 
reduce the remaining discrepancies by further in- 
creasing the statistics, however, would be a rather 
expensive enterprise. 

Let us therefore consider the cluster-diameter 
distribution which turned out to be a much 
cheaper solution of this problem. At first sight the 
constrained one-parameter fit to G in Fig. 2 
looks less perfect, since the data points are ran- 
domly scattered around the fit. The reason is that 
the correlations between the estimates at x and 
x + Ax are much smaller for G """ than for . 
This can be understood by noting that a cluster 
of diameter Xq contributes only to the one esti- 
mate of G (x) at x = xo, but to all estimates 
of g(°) (x) with x < Xo (recall the cluster estima- 
tor (5)). The correlation lengths resulting from 
unconstrained fits to G """ are collected in Ta- 
ble 1 . We see that the results are now in very good 
agreement with the theoretically expected values, 
with slight systematic deviations of only about 
1% — 2%. Contrary to g^(x) the fitted values 
are now overestimates. This tendency becomes 
obvious in Fig. 3 where we show the effective cor- 
relation lengths £f = l/ln[C(a;)/C(a; + 1)], with 
C = G or g(°\ Recall that neighboring val- 



ues of G are much less correlated than those 
of g(°K This explains the much larger error bars 
on the data from G . For q = 15 and 20 we 
observed a similar behaviour but, as is demon- 
strated in Fig. 3 for q = 20, here the £f of G d " m 
develop a much more pronounced plateau before 
the statistical errors increase and the data start to 
fluctuate around the theoretically expected value. 

Further investigations of the Ising and 3-state 
Potts model revealed, however, that the new ob- 
servable is not always advantageous. Our results 
for the Ising model in the disordered phase at 
/3 = 0.70340888 « 0.8/3 c are shown in Fig. 4. 
Here the £^ ff derived from G clearly over- 
shoot the exact value of ^ = 2.6202906 . . . before 
they slowly approach it from above. The Q s of 
g(°\ on the other hand, coincide with the exact 
value already for very small x, and a simple two- 
parameter fit of the form (6) with b = c = in 
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Figure 3. Effective correlation lengths for q = 10 
and 20 at (3 t in the disordered phase derived from 
the correlation functions as shown in Fig. 2. 
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Table 1 

Correlation length £,d{Pt) in the disordered phase extracted from the cluster-diameter distribution 
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300 x 300 






120 x 120 






80 x 80 


I min -I ml 




X min 


-x n 


Zd 


x min 
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40-130 


10.90(2) 
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-50 


4.297(4) 


13- 


-40 


2.766(3) 


64-130 


10.84(5) 


29- 


-50 


4.26(2) 


19 


-40 


2.744(7) 


88-130 


10.5(2) 


38- 


-50 


4.25(4) 


25- 


-40 


2.70(3) 


exact 


10.559519... 






4.180954. . . 






2.695502. . . 



the range x = 1...40 yields £ d = 2.62029(14), 
in excellent agreement with the exact value. The 
results for the 3-state model look very similar. 

Also for the three-dimensional g-state Potts 
model, which undergoes a first-order phase tran- 
sition already for q > 3, plots of the effective cor- 
relation lengths derived from G and g^ look 
qualitatively as in Fig. 4. Here we simulated the 
models with q = 3, 4, and 5 at the transition 
point p t in the disordered phase. Also in these 
cases we found that g^ gives much more reliable 
estimates of 

4. DISCUSSION 

Our numerical results clearly show that the 
cluster-diameter distribution is very well suited 
to extract the correlation length (,d(Pt) of two- 
dimensional g-state Potts models with relatively 
large values of q. While analyses of the stan- 
dard (projected) two-point function are plaqued 
by large systematic errors, with the new observ- 
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Figure 4. Effective correlation lengths for the 
Ising model at /3 = 0.70340888. 



able we succeeded to reproduce the theoretically 
expected values at a 1% — 2% level. 

For small values of q and also for the three- 
dimensional model at the first-order transition 
point, however, the standard correlation function 
gives much more reliable results. For reasons 
not well understood to date, the two quite dif- 
ferent correlators seem to behave complementary 
to each other. 
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